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Abstract 

The interplay between the collective dynamics of the quadrupole and octupole deformation de- 
gree of freedom is discussed in a series of Sm and Gd isotopes both at the mean field level and 
beyond, including parity symmetry restoration and configuration mixing. Physical properties like 
negative parity excitation energies, El and E3 transition probabilities are discussed and com- 
pared to experimental data. Other relevant intrinsic quantities like dipole moments, ground state 
quadrupole moments or correlation energies associated to symmetry restoration and configuration 
mixing are discussed. For the considered isotopes, the quadrupole-octupole coupling is found to be 
weak and most of the properties of negative parity states can be described in terms of the octupole 
degree of freedom alone. 



PACS numbers: 21.60.Jz, 27.70.+q, 27.80,+w 



I. INTRODUCTION. 



The nuclear mass region with proton number Z 60 and neutron number iV m 90 



ly, since it is a 
l). Particularly 



is receiving at present much attention, both experimental and theoretica 
region where nuclear structure collective effects of different nature overlap 
interesting in this context is the interplay between quadrupole transitional properties in 
N « 90 isotones and octupole deformation manifestations in nuclei with proton Z w 56 and 
neutron N pa 88 numbers. On one hand, isotones with N ^ 90 have been found as empirical 
realizations [2] of the critical point symmetry X(5), introduced to describe analytically 
the first order phase transition from spherical [U (5)] to well deformed [S77(3)] nuclei. Such 
critical point symmetries, have recently been studied within various microscopic approaches, 
either relativistic or non-relativistic (see, for example, and references therein). 

On the other hand, it is well known [1] that there is a tendency towards octupolarity 
around particular neutron/proton numbers, namely N/Z = 34, 56, 88 and 134. The emer- 
gence of octupolarity in these nuclear systems can be traced back to the structure of the 
corresponding single-particle spectra which exhibit maximum coupling between states of 
opposite parity, where the (N + 1, 1 + 3,j + 3) intruder orbitals interact with the (N,l,j) 
normal-parity states through the octupole component of the effective nuclear Hamiltonian. 
When the mixing is strong enough, the nucleus displays an octupole deformed ground state 
In particular, for nuclei with Z m 56 (N pa 88) the coupling between the proton 
(neutron) single-particle states /in/2 (213/2) and (/V/2) has been considered as mainly 
responsible for mean field ground state octupolarity. 

The search for signatures of stable octupole deformations in atomic nuclei has been ac- 
tively pursued d Uralg the last deeades flEJ.Asa m ain t ea tur e, oetupole defend evea-evea 
nuclei display particularly low-lying negative-parity 1~ states. In the case of stable octupole 
deformations, the + and 1~ states represent the members of parity doublets, giving rise to 
alternating-parity rotational bands with enhanced El transitions among them. These finger- 
prints of octupole deformations have already been found in the particular regions mentioned 
above, but especially in the rare-earth and actinide regions ilia]. 

For the sample of nuclei considered in the present study (i.e., 146_154 Sm and 148_156 Gd), 
experimental fingerprints have been obtained through the observation of octupole correla- 
tions at medium spins, as well as the crossing of the octupole and the ground state band, 



2 





FIG. 1: MFPESs computed with the Gogny-DlS EDF in panel a) and Gogny-DIM EDF in panel 
b) for the nucleus 150 Sm. Taking the lowest energy as a reference, the contour lines extend from 
0.25 MeV up to 4 MeV in steps of 0.25 MeV. Full, long dashed and short dashed contours are used 
successively to help identify contours more easily. Dotted lines correspond to contours starting at 
5 MeV and extending to 8 MeV in steps of 1 MeV. 



jointing to the fact that reflection symmetric and asymmetric structures coexist in Sm 



9| and 148 Sm [h]]. A recent study 11] has analyzed the lowest four negative-parity bands 
in 152 Sm and has found an emerging pattern of repeating excitations, built on the 0^ level 
and similar to that of the ground state, suggesting a complex shape coexistence in 152 Sm. 
The experimental findings 9Mll| mentioned above, already suggest that it is timely and 



necessary to carry out systematic studies of the quadrupole-octupole interp 
other regions of the nuclear chart, starting from modern (global) relativistic 



non-relativistic 



13 



ay 


in 


12. 


13 



131 ] and/or 



16j nuclear Energy Density Functionals (EDFs), with reasonable predic- 
tive power all over the nuclear chart. 

Let us remark that the microscopic study of the dynamical (i.e., beyond mean field) 
quadrupole-octupole coupling in the considered Sm and Gd isotopes is also required to 
better understand the extent to which a picture of independent quadrupole and octupole 
excitations persists or breaks down for nuclei with neutron number ^ ^ 88. This, together 
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TABLE I: Proton E P (Z) (MeV) and neutron E p (N) (MeV) pairing energies, dipole D (efm) 
moment, quadrupole Q2o(b) an d octupole Q3o(b 3 ^ 2 ) moments at the minima of the MFPESs for 
the isotopes 146 ~ 154 Sm. Results are given for both Gogny-DlS and Gogny-DIM EDFs. 



Nucleus 


E p (Z) E p (N) 


Dq Q20 Q30 


E p (Z) E p (N) 


Dq Q20 Q30 






D1S 




DIM 


146 Sm 


-14.64 


-3.58 


0.00 1.20 0.00 


-15.67 


-5.20 


0.00 0.60 0.00 


148 Sm 


-12.41 


-3.34 


0.18 3.00 1.25 


-13.92 


-4.52 


0.14 3.00 0.75 


150 Sm 


-10.98 


-1.53 


0.41 4.80 1.50 


-11.80 


-2.99 


0.35 4.80 1.25 


152 Sm 


-6.58 


-5.57 


0.00 7.20 0.00 


-8.24 


-6.18 


0.00 6.60 0.00 


154 Sm 


-5.91 


-3.27 


0.00 7.80 0.00 


-6.38 


-4.63 


0.00 7.80 0.00 



with the available experimental fingerprints |9Hll| for octupolarity in the region, is one of 
the main reasons driving our choice of the nuclei 146_154 Sm and 148_156 Gd as a representative 
sample to test the performance of the different approximations and EDFs considered in the 
present study. 

From a theoretical perspective, many different models have been used to describe octupole 
correlations in atomic nuclei. For a detailed survey the reader is referred, for example, to Ref. 
[l| . Calculations based on the shell-correction approach with folded Yukawa deformed poten- 
tials [17, 181, as well as calculations based on Woods-Saxon potentials with various models 

nn 

for the microscopic and macroscopic terms |19l. 1201]. predicted a significant stabilization of 
octupole deformation effects in various nuclear mass regions. Pioneer Skyrme-HF+BCS 
calculations including the octupole constraint and restoring parity symmetry were carried 
out in Ref. {2]]]. Subsequent calculations in Ref. {22]] included both quadrupole and oc- 
tupole constraints at the same time but at the mean field level only. On the other hand, 
microscopic studies of octupole correlations with Skyrme and Gogny EDFs, both at the 



mean field leve 
(see, Refs 



and beyond with different levels of complexity, have already been reported 
23f|33| and references therein) for several regions of the nuclear chart. The- 
oretical studies in the Sm region include mean field based calculations with the collective 
hamiltonian and the Gogny force [3^, the IBM study with spdf bosons of Ref. 34] or the 



collective models using a coherent coupling between quadrupole and octupole modes 



35] 



and new parametrizations of the quadrupole and octupole modes [36|. Non-axial pear-like 
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shapes in this region were considered, for example, in Refs. [371]. Additionally, the isotopes 
i46-i56g m have been investigated very recently within the constrained reflection-asymmetric 
relativistic mean field (RMF) approach 38[ based on the parametrization PK1 [39j for the 
RMF Lagrangian together with a constant gap BCS approximation for pairing correlations. 

In the present work, we investigate the interplay between octupole and quadrupole de- 
grees of freedom in the sample of nuclei 146_154 Sm and 148_156 Gd. We use three different levels 
of approximation. First, the constrained (reflection- asymmetric) Hartree-Fock-Bogoliubov 
(HFB) framework is used as starting point providing energy contour plots in terms of the 
(axially symmetric) quadrupole Q 20 = (^l^ol^) an d octupole Q 30 = (<&\Q 30 \$) moments 
(where |$) is the corresponding HFB intrinsic wave function). Within this mean field frame- 
work we pay attention to the shape changes in the considered nuclei and their relation with 
the underlying single-particle spectrum 



32 



40]. 



As will be discussed later on, the (Q20, Q30) mean field potential energy surfaces (MF- 
PES) obtained for the nuclei 146_154 Sm and 148_156 Gd are, in most of the cases, very soft 
along the Q30 direction indicating that the (static) mean field picture is not enough and that 
a (dynamical) beyond mean field treatment is required. Therefore, both the minimization 



of the energies obtained after parity projection of the intrinsic states |21l . 1301 . |33| as well as 
quadrupole-octupole configuration mixing calculations in the spirit of the Generator Coor- 
dinate Method (GCM) |4l|, are subsequently carried out. The analysis of the two sets of 
results allows to disentangle the role played in the dynamics of the considered nuclei by the 
restoration of the broken reflection symmetry and the fluctuations in the (Q2o,Q3o) collective 
coordinates. Similar calculations with the Skyrme functional where carried out in Ref 26 1 
for a lead isotope. 

To the best of our knowledge, the hierarchy of approximations (i.e., reflection-asymmetric 
HFB, parity projection and (Q20, Qso)-GCM) considered in the present work belong, at least 
for the case of the Gogny-EDF, to the class of unique and state-of-the-art tools for the mi- 
croscopic description of quadrupole-octupole correlations in atomic nuclei. Let us also stress 
that, the two-dimensional GCM (2D-GCM) framework used in the present study repre- 



sents an extension of the treatment of octupolarity reported in Refs. |32l . |33| . where a 
one- dimensional collective hamiltonian based on several approximations and parameters ex- 
tracted from (^so-constrained HFB calculations, was considered. Here, on the other hand, 
the octupole and quadrupole degrees of freedom are explored simultaneously and the kernels 
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TABLE II: The same as Table I but for the isotopes 



148-156 



Gd. 





E p (Z) E p (N) 


-Do Q20 Q30 


E p (Z) E p (N) 


Dq Q20 Q30 






D1S 




DIM 


148 Gd 


-15.22 


-4.27 


0.00 0.66 0.00 


-16.11 


-5.35 


0.00 0.00 0.00 


150 Gd 


-14.22 


-3.63 


0.19 3.60 0.75 


-15.43 


-5.03 


0.05 3.00 0.25 


152 Gd 


-12.69 


-3.02 


0.27 4.80 1.00 


-13.18 


-4.76 


0.15 4.80 0.50 


154 Gd 


-7.63 


-6.26 


0.00 7.20 0.00 


-9.25 


-6.88 


0.00 6.60 0.00 


156 Gd 


-7.18 


-4.86 


0.00 7.80 0.00 


-7.66 


-6.20 


0.00 7.80 0.00 



involved in the solution of the corresponding Hill- Wheeler equation [41[ are computed with- 
out assuming a gaussian behavior of the norm overlap neither a (second order) expansion 
over the non-locality of the hamiltonian kernel. Therefore, the present study for the selected 
set of Sm and Gd nuclei, to the best of our knowledge the first of this kind for the case of the 
Gogny-EDF, may also be regarded as a proof of principle concerning the feasibility of the 
calculations to be discussed later on. Pioneer calculations along the same lines considered 
in the present study, but based on the Skyrme-EDF, have been carried out in Ref. 



26 



In addition to the standard Gogny-DlS 



42]. 



15] parametrization, which is taken reference, 



the DIM parametrization 43j will also be considered. The functional Gogny-DlS has a long 
standing tradition and it has been able to describe many low-energy experimental data all 



over the nuclear chart with reasonab 
beyond (see, for example, Refs. [la. 



e predicti ve p ower both at the mean field level and 



27 



31 



44 



541 ] and references therein). On the other 



hand, the DIM parametrization [43j that was tailored to provide a better description of 
masses is now proving its merits in nuclear structure studies not only in even-even nuclei 



40 



43 



(EFA) 



544 
54 



58 
56 



, but also in odd nuclei in the framework of the Equal Filling Approximation 
58| . In this paper the results of both D1S and DIM are compared to verify 



the robustness of our predictions with respect to the particular version of the interaction 
and to test the performance of DIM in the present context of quadrupole-octupole coupling. 

The paper is organized as follows. In Sees. [IJ],[TIIJand[rV]we briefly describe the theoretical 
formalisms used in the present work and subsequently the results obtained with them. Mean 
field calculations will be discussed in Sec. (TTJ Parity projection and configuration mixing 
results will be presented in Sees. [TTTJand[IVl respectively. In particular, in Sec. [TV] especial 
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FIG. 2: In panel a) the reflection symmetric (i.e., Q 30 =0) MFPECs for 146 ~ 154 Sm and in panel b) 
for 148 - 156 Gd are plotted as functions of the axially symmetric quadrupole moment C^o- Results 
for both Gogny-DlS with full line and Gogny-DIM with dashed line are given. In each panel the 
energies are referred to the D1S ground state energy of the heavier isotope. 

attention will be paid to beyond mean field properties in the considered nuclei -dynamical 
octupole and dipole moments, correlation energies, reduced transition probabilities B(E1) 
and B(E3) as well as energy splittings- and their comparison with available experimental 
data. Finally, Sec. |V] is devoted to the concluding remarks and work perspectives. 



II. MEAN FIELD SYSTEMATICS. 

The aim of the present work is the study of the quadrupole-octupole dynamics in se- 
lected Sm and Gd isotopes with neutron number 84 < N < 92. Three different levels of 
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FIG. 3: (Color online) Single particle energies (see text for details) in 150 Sm are plotted as a 
function of the quadrupole moment Q20 (in barn) up to the value corresponding to the ground 
state minimum. From there on, the plot continues with the representation of the SPEs as a function 
of the octupole moment Q30. In the part of the plot where the SPEs are plotted versus Q20, full 
(dashed) curves stand for positive (negative) parity levels. The thick (red) dashed line in each plot 
represents the chemical potential. In the part of the plot where the SPEs are plotted versus Q30 
some levels around the Fermi level are labeled with twice their J z value. Panels a) and b) (c) and 
d)) correspond to results obtained with D1S (DIM) EDFs. Panels a) and c) (b) and d)) correspond 
to protons (neutrons). 



approximation are considered: the HFB method with constraints in the relevant degrees 
of freedom, parity projection (with minimization of the energy after projection) and the 
Generator Coordinate Method (GCM) with Q 2 o and Q 30 as collective coordinates. For a 
detailed survey on the three techniques the reader is referred to Ref. {4]]]. The Gogny EDF 
is used consistently in the three methods both with the D1S and DIM parametrizations. 

First, (Q20, Q3o)-constrained HFB calculations are performed for the nuclei 146_154 Sm and 
148_156 Gd to obtain a set of states |$(Q)) labeled by their corresponding multipole moments 
Q = (Q2o,Q3o)- The K — quadrupole Q20 and octupole Q30 moments are given by the 
average values 

Q2o = (^ 2 - 1 -(x 2 + y 2 )\^} } (1) 



Q, = ^-l( x ^ + y A z \^). (2) 



and 

Axial and time reversal are self-consistent symmetries in the mean field calculations. As a 
consequence of the axial symmetry imposed on the HFB wave functions |$), the mean values 
of the multipole operators Q211 and with /x 7^ are zero by construction. Aside from the 
constraints on the quadrupole and octupole moments, a constraint on the center of mass 
operator is used to place it at the origin of coordinates in order to prevent spurious effects 
associated to center of mass motion. The HFB quasiparticle operators 41 j have 

been expanded in an axially symmetric harmonic oscillator (HO) basis (cj\q) containing 13 
major shells as to grant convergence for all the observable quantities. For the solution of 
the HFB equation, an approximate second order gradient method [59(] is used. 

The MFPES have been computed in a grid with Q20 in the range from -30b to 30b in 
steps of 0.6 b and the octupole moment Q30 in the range from b 3 / 2 to 3.75 b 3 / 2 in steps 
of 0.25 b 3 / 2 . Negative values of the octupole moment are not computed explicitly as the 
corresponding wave function can be obtained from the positive Q™ one by applying the 
parity operator. As the Gogny EDF is invariant under parity (see 44j, |45( for a discussion 
of the meaning of symmetry invariance for density dependent "forces") the energy has the 
property Ehfb(Q2o,Q3o) = Ehfb(Q2o, —Q30) and therefore is an even function of the oc- 
tupole moment. For this reason, in the graphical representation of the PES only positive 
values of Q 30 are considered. 

The MFPESs obtained for the nucleus 150 Sm, with the parametrizations D1S and DIM 
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FIG. 4: Positive 7r=+l (upper panels) and negative tt=-1 (lower panels) parity-projected potential 
energy surfaces (PPPES) computed with the Gogny-DlS (left panel) and Gogny-DIM (right panel) 
EDFs for the nucleus 150 Sm. See caption of Fig [T] for the contour line patterns. 

of the Gogny-EDF, are shown in Fig. [I] as an illustrative example of our mean field results. 
For the sake of presentation, quadrupole and octupole moments have been constrained in 
the plots to the ranges -10 b < Q 20 < 20 b and b 3 / 2 < Q 3 o < 3.75 b 3 / 2 , respectively. 
The similitude between the D1S and DIM results both in the Q20 and Q30 directions is 
remarkable. In previous calculations in other regions and looking at different physical effects 



54 



56N58| we have already noticed the same similitude between D1S and DIM results. 
Focusing on the MFPES, the absolute minimum is located in the prolate side at a finite value 
of the octupole moment. The minimum is very shallow along the Q30 direction. Another 
minimum is observed in the oblate side, but this time centered at Q30 = 0. For the other 
nuclei considered the energies look similar and therefore they are not shown. The most 
relevant mean field quantities for the ground states are summarized in Tables I and II. In 
order to better understand the quadrupole deformation properties of the studied nuclei, 
the reflection symmetric (i.e., Qso=0) mean field potential energy curves (MFPECs) are 
depicted for all the considered nuclei in Fig. [2j A transition from weakly deformed ground 
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FIG. 5: Mean field (full line), positive (dashed) and negative (dotted) parity-projected energies as 
a function of the octupole moment Q30 for selfconsistent Q20 values for the nucleus 150 Sm. Energies 
are referred to the D1S HFB ground state energy. 



states in the N=84 nuclei 146 Sm and 148 Gd to well (quadrupole) deformed ground states in 
i52,i54g m anc j 154,1560^ (p ro ] a te moments 6.6 b < Q20 < 7.8 b) is observed. In most of the 
isotopes except the lightest ones an additional minimum is observed in the oblate side. This 
minimum may become a saddle point (see, [6] for examples) once the 7 degree of freedom is 
considered. Nevertheless, the simultaneous consideration of triaxial quadupole and octupole 
moments lies outside of the scope of the present study. Investigation along these lines is in 
progress and will be reported elsewhere. 
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From Tables I and II, we observe the onset of an octupole deformed regime at the N=88 
nuclei 150 Sm and 152 Gd. These nuclei mark the borders of another shape transition from 
octupole deformed ground states in 148 Sm and 150 Gd to quadrupole deformed and reflection 
symmetric ground states in 152 Sm and 154 Gd. Consistent with the breakdown of the left- 
right symmetry in their ground states, the 148 > 150 Sm and 150 ' 152 Gd isotopes exhibit a non 
zero (static) dipole moment Dq. It is computed as the ground state average value of the 
dipole operator 

d --z -?-z m 

along the symmetry z-axis. The values of D Q tend to be smaller for DIM than for D1S. 
This is not surprising due to the delicate balance between single particle orbital properties 



that enter in the definition of the dipole moment 29|]. Another quantity of interest is the 
mean field octupole correlation energy E^ r = E g HFB Q 3Q=0 — E 9 H S FB corresponding to the 
energy gain by allowing octupole deformation. For example, the values obtained for 150 Sm 
and 152 Gd are 204 and 43 KeV (105 and 6 KeV) for the functional D1S (DIM), respectively 
These very low values are a clear indication of the softness of the octupole minima in those 
nuclei. As the minima are also soft along the Q20 direction both the quadrupole and octupole 



degrees o 
problem 



: freedom have to be considered at the same time in a dynamical treatment of the 



sHii|. 



In Tables I and II, the proton E p (Z) and neutron E p (N) pairing energies are also listed. 
They are computed in the usual way as E p {r) = -l/2Tr(A( t)k*(t)) in terms of the pairing 
field A and the pairing tensor k for each isospin r = Z, N. Moving along isotopic chains, 
the smallest neutron pairing energy corresponds to the N=88 nuclei 150 Sm and 152 Gd, which 
are precisely the ones providing the largest values of the mean field octupole correlation 
energy E^ r . The significant lowering of the neutron pairing energies in these nuclei is a 
consequence of the low level density typical of deformed (quadrupole or octupole) minima, 
the Jahn- Teller effect. On the other hand, proton pairing energies tend to decrease as a 
function of the neutron number. In general, the proton and neutron pairing energies for 
the two Gogny-EDFs considered follow the same trend, the only relevant difference being in 
their absolute values that tend to be slightly larger for DIM. 

Before concluding this section, we turn our attention to single-particle properties. The 
appearance of quadrupole and/or octupole deformation effects is strongl y l inked to the 



position of the Fermi energy in the single-particle spectrum 
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FIG. 6: Collective wave functions squared (\G a (Q 20, Q30) | 2 ) for the ground state (left panels) and 
the lowest negative parity state (right panels) for the nuclei 46 ' 150,154 Sm. ^he contour lines (a 
sucession of full, long dashed and short dashed) start at 90 % of the maximum value up to 10 % 
of it. The two dotted line contour correspond to the tail of the amplitude (5 % and 1% of the 
maximum value). 

the evolution of the the single-particle energies (SPEs) for both protons and neutrons with 
deformation is an interesting piece of information. In HFB calculations the concept of single 
particle energy is assigned to the eigenvalues of the Routhian h = t + Y — Aq 20 Q 2 o — ^q-mQso, 
with t being the kinetic energy operator, and T the Hartree-Fock field. The term Aq 20 <520 + 
Aq 30 Q3o contains the Lagrange multipliers used to enforce the corresponding quadrupole and 
octupole constraints. 

Proton and neutron SPEs for the nucleus 150 Sm, computed with both the Gogny-DlS 
and Gogny-DIM EDFs are presented in Fig. [31 The SPEs are plotted first as functions of 
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TABLE III: Dynamical quadrupole Q$ , Q { 20 > (b) and octupole Q { 3 q> , Q { 30 > (b 3 / 2 ) moments 
corresponding to the first positive and negative parity 2D-GCM states in the isotopes 146_154 Sm. 
Results are given for both Gogny-DlS and Gogny-DIM EDFs. 

Nucleus Q^q O30 ^20^ ^?30^ Qso^ Qw^ O30 

DIS DIM 



146 Sm 


0.63 


0.43 


1.37 


1.33 


0.45 


0.39 


1.39 


1.29 


148g m 


2.28 


0.54 


2.01 


1.77 


2.93 


0.52 


2.10 


1.65 


150 Sm 


2.94 


0.60 


2.63 


1.83 


3.09 


0.56 


2.68 


1.81 


152 Sm 


5.48 


0.51 


3.81 


1.72 


5.28 


0.50 


3.63 


1.74 


i54 Sm 


6.15 


0.50 


4.21 


1.63 


6.15 


0.49 


4.33 


1.58 



TABLE IV: The same as Table III but for the isotopes 148 ~ 156 Gd. 



Nucleus 


^20 


o {+) 


^20 


^30 


^20 


{+) 
V30 


V20 


V30 






DIS 






DIM 




148 Gd 


0.23 


0.44 


1.05 


1.35 


0.12 


0.41 


0.97 


1.29 


150 Gd 


2.46 


0.52 


1.78 


1.74 


2.53 


0.47 


1.57 


1.65 


152 Gd 


3.47 


0.57 


2.66 


1.79 


3.50 


0.55 


2.73 


1.73 


154 Gd 


5.72 


0.50 


3.75 


1.74 


5.50 


0.49 


3.31 


1.68 


156 Gd 


6.51 


0.49 


4.47 


1.59 


6.40 


0.48 


4.56 


1.61 



the quadrupole moment Q20 up to the value corresponding to the ground state minimum 
obtained with the Q30 — constraint. From there on, the plot continues with the represen- 
tation of the SPEs as a function of the octupole moment Q^. The given SPEs as a function 
of the octupole moment have the self-consistently determined quadrupole moment which, in 
the present case, does not depart significantly from the ground state value at Qzo = 0. 

The first significant conclusion drawn from Fig. [3] is that the DIS and DIM SPE plots 
look rather similar near the Fermi level (thick red dashed line): both the ordering of the 
levels at sphericity and their behavior with Q 2 o and Q30 are rather similar. For this reason 
we will from now on focus only on the DIM SPEs. For protons, the positive parity 
orbital strongly interacts with the negative parity /in/2 one by means of the I = 3 octupole 
component of the interaction. The position of the proton's Fermi level in the considered 
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nucleus is located in the center of a small gap in the single particle spectrum that favors 
octupole deformation (Jahn- Teller effect 621]) . In the neutron's spectrum a fairly large 
gap near the Fermi level also opens up when the octupole moment is switched on. The 
neighboring levels come from the negative parity fj/2 orbital and the positive parity 113/2 
intruder orbital. It is also worth mentioning the occurrence of "quasi-j" orbitals in the 
neutron spectrum for the Q30 values corresponding to the minimum at around 2 b 3//2 . A j = 
7/2 is formed at an energy of around —4 MeV; one with j = 5/2 is located at around —6 MeV 
and finally another one with j =3/2 shows up at an energy of —8 MeV. The same grouping 
of levels can also be observed in the SPEs for protons at similar values of the octupole 
moment. These quasi- j orbitals are the consequence of the relationship between classical 
closed periodic orbits for specific octupole deformed shapes and the corresponding quantum 
orbitals that have to show an integer ratio between the radial and angular frequencies (see 



6l| , Vol II, page 587 for a general discussion and also [29j for specific examples in rare earth 
nuclei). 



III. PARITY PROJECTION. 



Although the HFB framework discussed in the previous section is a valuable starting 
point, it produces MFPESs with very soft minima along the Qso-direction in the nuclei con- 
sidered. This suggests the important role played by both types of dynamical correlations: 
the one associated with symmetry restoration and the other to configuration mixing. Sym- 
metry restoration is considered in this section while configuration mixing will be presented 
in the next section. 

There are two spatial symmetries broken in the present calculations. One is rotational 
symmetry with the quadrupole moment as relevant parameter and the other is reflection 
symmetry (parity) with the octupole moment as relevant quantity. From the discussion of 
the mean field results it is clear that the softest mode is the octupole moment and therefore 
the most relevant symmetry to be restored is parity. Obviously, it would be desirable to 
restore also the rotational symmetry as well as particle number. This combined symmetry 
restoration is feasible but, when combined with the configuration mixing of next section, 
becomes a very demanding computational task not considered in this paper. 

The quantum interference typical of the GCM framework could be directly used to restore 
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FIG. 7: (Color online) The dynamical dipole moments provided by the GCM calculations for the 
nuclei 146_154 Sm and 148_156 Gd are shown, as functions of the neutron number, in panels a-b). 
Experimental dipole moments are taken from Ref. [II] . Results are shown for both Gogny-DlS and 
Gogny-DIM EDFs. 

the parity symmetry by choosing appropriate weights for the configurations with mutipole 
moments (Q2o,Qso) and (Q20, —Q30) [33f . However, in order to disentangle the relative 
contribution of the parity restoration correlations as compared with the ones of the GCM 
configuration mixing, we have carried out explicit parity projection calculations. 

To restore parity symmetry {21} [3^] we build positive (n = +1) and negative ( 71 = —1) 
parity-projected states |$ 7r (Q2o, Q30)) = 'P n \^(Q2o, Q30)) by applying the parity projector 
V 77 to the intrinsic configuration. The parity projector is a linear combination of the identity 
and the parity operator II given by 



V 



(4) 



The projected energies, used to construct parity-projected potential energy surfaces (to be 
called PPPES in what follows), are labeled with the multipole moments Q = (Q20, Q30) of 
the intrinsic state and read 63] 

E(Q)= (HQ)\H[ P (mm)) ($(Q)|W)]II|$(Q)) 

^ <$(Q)|$(Q))+7r($(Q)|n|$(Q)> ($(Q)|$(Q)) + 7r($(Q)|n|$(Q)> 



The parity projected mean value of proton and neutron number, ^(q|^-|^q)^ an d 
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FIG. 8: Predicted energy splittings between the lowest lying ir=+l and tt=-1 2D-GCM states in 
146-I54g m anc j 148-1560^ are CO mpared with the experimental + — 1~ and + — 3~ splittings 



Results are shown for the D1S parametrization of the Gogny force as DIM ones are rather similar. 



^HQ^^l^iQ)) usua ^y differ from the nucleus' proton Z and neutron N numbers. To 
correct the energy for this deviation we have replaced H by H — \ z (z — Z^j —X n (n — N$_ , 
where Xz and Xn are chemical potentials for protons and neutrons, respectively 3fj, |64], 165 ]. 

In the case of the Gogny-EDF, as well as for Skyrme-like EDFs, the definite expression 
for the projected energy (JSJ) depends, on the prescription used for the density dependent 
part of the functional. In this work, we resort to the so called mixed density prescription 
that amounts to consider the standard intrinsic density 

($(Q)|p(f)|$(Q)> 



p(r) 



and the density 



<$(Q)|$(Q)> ' 
($(Q)|p(r)n|$(Q)) 



(6) 



(7) 



($(Q)|IT|$(Q)> 

in the evaluation of the first and second terms in Eq. (J5]), respectively. The mixed density 
prescription has been widely and successfully used in the context of projection and/or con- 
figuration mixing techniques (see, for example, 13|, |44], |45|, |65M68| and references therein). 



In fact, this is the only prescription that guarantees various consistency requirements within 
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the EDF framework 



put into evidence recently 



69| . Even though this prescription has some drawbacks, as 



69h71|. the use of other prescriptions, like the one based on the 



projected density, are pathologically ill defined when applied to the restoration of spatial 
symmetries j^ . 

As an illustrative example of PPPES, we show in Fig. S]the results for the nucleus 150 Sm 
obtained with both the D1S and DIM parametrizations of the Gogny force. Along the 
Q30=0 axis, the projection onto positive parity ir=+l is unnecessary since the corresponding 
(quadrupole deformed) intrinsic configurations are already parity eigenstates with eigenvalue 
7r = +1. For the same reason, the negative parity n = — 1 projected wave function only 
makes sense along the Qso=0 axis when a limiting procedure is considered. The evaluation of 
physical quantities in this case is subject to numerical inaccuracies consequence of evaluating 
the ratio of two small quantities (the denominator is the norm of the projected negative parity 
state that is zero in this case) and alternative expressions, obtained by considering explicitly 



the Q3o= limit [30[, are required for a sound numerical evaluation of those quantities. 



Note however (see, Fig. |5]) that the negative parity projected energy increases rapidly while 
approaching the Q3o=0 configuration and therefore it does not play a significant role in the 
subsequent discussion of the corresponding PPPESs. As a consequence, we have omitted 
this quantity along the Q3o=0 axis. 

As in the mean field case, the results with D1S and DIM show a striking similarity and 
therefore only the D1S results will be discussed. The comparison between the MFPESs 
in Fig. [1] and the PPPESs in Fig. HJ clearly illustrates the topological changes induced 
by the restoration of the reflection symmetry. In general, the quadrupole moments Q20 
corresponding to the absolute minima of the PPPESs, remain quite close to the ones obtained 
at the HFB level (see, Tables I and II) increasing their values as more neutrons are added 
for each of the Sm and Gd chains. On the other hand, the situation is quite different along 
the Q30 direction. To obtain a more quantitative understanding of the evolution of the 
PPPESs, we have plotted in Fig. [5] the parity-projected energy curves for selfconsistent Q20 
values, as a function of the octupole moment Q30 for the nucleus 150 Sm. The corresponding 
HFB energy curves are also included for comparison. For 150 Sm, and all the other nuclei 
considered in the present study, the negative parity curves always show a well developed 
minimum at Q30 values in the range 1.50-1.75 b 3 / 2 . On the other hand, and regardless of 
the particular version of the Gogny-EDF employed, the 7r = +1 curves always display a 
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characteristic pocket 



30] with a minimum at Q30 = 0.50-0.75 b 3 / 2 . In the spirit of the 
Variation after Projection procedure, the configuration yielding the minimum of the positive 
(negative) parity projected energy as a function of Q 2 o and Q 30 is to be associated with the 
positive (negative) parity state. As a consequence of this "minimization after projection" 
the intrinsic states for each parity have different deformations. The positive parity ground 
state gains an amount of energy EFJrr ^ T °^ given by 

z-.par proj _ rpg.s. j^g.s. , & \ 

Morr - ^HFB ~ K=+l (°J 

where, E%'1' +1 corresponds to the absolute minima of the positive parity PPPESs and E 9 H S FB 
to the HFB ground state energies, i.e., the absolute minima of the MFPESs. Regardless of 
the Gogny-EDF employed, they are always smaller than 900 keV in each of the considered 
nuclei. This correlation energy has to be compared to the correlation energy gained by 
configuration mixing (see also Fig. [TOl below) . 

IV. GENERATOR COORDINATE METHOD. 

According to the discussions in previous sections, it can be concluded that not only the 
plain HFB results of Sec. [TTJ but even the parity projection ones, may not be sufficient 
to decide whether, as suggested in Ref. [381], there exists a transition to an octupole de- 
formed regime in the considered nuclei in addition to the transitional behavior along the 
(520-direction {5, 6|. Within this context, {Q201 <?3o)-GCM calculations are needed in order 
to verify the stability of the quadrupole and/or octupole deformation effects encountered in 
both the MFPESs and the PPPESs for the considered Sm and Gd nuclei. One should also 
keep in mind that in the framework of such a dynamical 2D-GCM treatment, not only the 
mean field energy surface but also the underlying collective inertia plays a role. 

The superposition of HFB states 

TO = /dQ/;(Q)|$(Q)> (9) 
is used to define the GCM wave functions l^)- I n the integration domain both positive and 



negative octupole moments Q 30 are inc 



of the Hill- Wheeler (HW) equation [41] 



uded. The GCM amplitudes f£(Q) are the solutions 



J dQ' (H(Q, Q') - E%tf(Q, Q')) /;(Q') = 0. (10) 
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FIG. 9: Theoretical and experimental transition rates B(E1, 1~ — > + ) (panels a and b) and 
B(E3,3~ — > + ) (panels c and d) for the nuclei 146 ~ 154 Sm and 148 ~ 156 Gd. Results are shown for 

► + ) rates are extracted 



the Gogny-EDFs D1S and DIM. Experimental results for B(E1,1 
from Ref. [1] while the experimental B(E3,3~ — > + ) rates are taken from Ref. 



73]. 
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The GCM hamiltonian "H(Q, Q ) and norm jV(Q, Q ) kernels are given by 

H(Q,Q') = (HQ)\H[P GCM (mm')), 

AT(Q,Q') = ($(Q)|$(Q')> (11) 

where in the evaluation of %(Q, Q ) the mixed density prescription is used 

($(Q)|p(f)|$(Q')> 

p (r)_ mQ)\m')) ■ { } 



As in the parity projection case the Hamiltonian kernel K(Q, Q ] 



30. 



is also supplemented 



64 



65] 



with first order corrections in both proton and neutron numbers 

The solution of the HW equation ( fTOl) provides the energies corresponding to the 
ground (a = 1) and excited (a = 2,3,...) states. The parity of each of these states is 
given by the behavior of /^(Q) under the Q30 —> — Q30 exchange. This is a consequence 



of the invariance under reflection symmetry of the GCM Hamiltonian kerne 



on the solution of Eq. ([TO]) , the reader is referred, for example, to Refs. 



41 



s. 



44 



or details 



65] . Since 



the |$(Q)) basis states are not orthogonal, the functions /J(Q) of Eq. §§§ can not be 
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44j) 



interpreted as probability amplitudes. One then introduces (see, for example, Refs. 
the collective wave functions 

GKQ) = J dQ'jV^Q, Q')/;(Q'), (13) 

which are orthogonal and therefore their modulus squared |G^(Q)| 2 has the meaning of a 
probability amplitude. It is easy to show that the parity of the collective wave functions 
G a (Q2o, Q30) under the exchange Q30 — > — Q30 corresponds to the spatial parity operation 
in the correlated wave functions built up from them. The inclusion of octupole correlations 
immediately restores the reflection symmetry spontaneously broken at the mean field level 
and grants the use of a parity label it for the GCM quantities. 

The collective wave functions of Eq. fflBl can be used to express overlaps of operators 
between GCM wave functions in a more convenient way 

(*:|6|^:> = / dQdQ'G? (Q)O(Q, Q')G^(Q'), (14) 

with the kernels 

0(Q,Q') = J dQ"dQ"'Af-^Q-Q")(Q"\0\Q"')Af-HQ"'-M') (15) 
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FIG. 10: (Color online) The 2D-GCM correlation energies for the nuclei 146_154 Sm (panel a) and 
i48-i56Qcl (panel b) are shown, as functions of the neutron number. The correlation energies 
stemming from the restoration of reflection symmetry are also included. Results are shown for 
both Gogny-DlS and Gogny-DIM EDFs. For more details, see main text. 

given in terms of the operational square root of the overlap kernel that is defined by the 
property 

A/-(Q; Q') = J dQ"Aft (Q; Q")A/^ (Q"; Q'). (16) 

The solution of Eq. ffTTJ]) allows the calculation of physical observables like the energy 
splitting between positive and negative parity states as well as B(E1) and B(E3) transition 
probabilities. In the present study time reversal symmetry is preserved and therefore only 
excited states with an average angular momentum zero can be accounted for. Genuine 1~ 
and 3~ states, on the other hand, will require to consider cranking HFB states 3l|, |41[, a 



calculation which is out of the scope of the present work. We assume here that the cranking 
rotational energy of the 1~ and 3~ states is much smaller than the excitation energy of 
the negative-parity bandhead and therefore it can be neglected. For the reduced transition 
probabilities B(E1, 1~ — > + ) and B(E3,3~ — > + ) the rotational model approximation for 
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K=0 bands has been used 

e 2 



2 



B(EX, \- 0+) = — (*r~ l \Ox\n=i l ) ' (17) 

where a corresponds to the first GCM excited state of negative parity. The electromagnetic 
transition operators (D\ and 0% represent the dipole moment operator of Eq. ([3]) and the 
proton component Q 30 prot °f the octupole operator, respectively. The evaluation of the 
overlap is carried out using Eq. (fT4j) . 

In Fig. E]the collective probability amplitude |G£(<520> Q30) | 2 of Eq. (1T3T) , obtained from 
the solution of the HW equation ffTUl) are plotted. As a typical example, results for the 
I46,i50,i54g m j so t p e g anc i Gogny-DlS EDF are presented. For other nuclei and Gogny 
parametrizations, the results look very similar. The left panels in Fig. [6] correspond to the 
ground state wave functions (i.e., a = 1 and 7r=+l) while the right panels correspond to 
the lowest-lying tt=-1 states o = 3 for 146 Sm and a = 2 for the others. 

The ground state collective probability amplitude \G*zf 1 (Q2o, Q30) | 2 reach a global max- 
imum at Q30 = pointing to the octupole-soft character of the ground states in 146_154 Sm. 
The spreading along the octupole direction is large for 150 > 154 Sm indicating octupole softness 
in these nuclei. For the negative parity collective wave functions the maximum is always 
located at a non zero value of Q^o as could be anticipated from the parity projection re- 
sults. For 150 ' 154 Sm the wave function spreads out farther along Q 30 than in previous cases 
in agreement with the octupole softness of their ground states. 

To have a more quantitative characterization of the collective wave functions we have 
computed mean values of relevant operators (see Eq. (fl4l) ). The first is the average of the 
quadrupole moment defined as 

(Q 20 )Z=(VZ\Q W \*Z). (18) 

For negative parity operators like the octupole or the dipole moment the above averages are 
zero by construction and therefore a meaningful averaged quantity has to be defined as 

01 = A f , dQdQ'Gl*(Q)0(Q, Q')GJ(Q') (19) 

where a restriction to positive values of the octupole moment has been made. The average 
quadrupole and octupole moments for the ground state (a = 1) are listed in 
Tables III and IV. The Q20 moments follow a trend similar to the one found within the 
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FIG. 11: (Color online) The sum (blue curve named 1D-Q2+1D-Q3) of the correlation energies 
Ecorr® 3 [Eq.©] (green curve named 1D-Q3) and E^~ Q2 [Eq.flMD] (red curve named 1D-Q2) is 
compared with the correlation energy E^®~ GCM [Eq. (j20p ] (black curve named 2D-GCM) provided 
by the full 2D-GCM calculations in 146_154 Sm and 148 ~ 156 Gd. Results are shown for both Gogny- 
D1S and Gogny-DIM EDFs. For more details, see main text. 

HFB approximation increasing their values as more neutrons are added in a given isotopic 
chain. On the other hand, the isotopic trend predicted for is quite different to the one 
predicted at the mean field level. As discussed in Sec. Ull at the Gogny-HFB level only the 
N=86 and 88 isotones 148 > 150 Sm and 150 ' 152 Gd display non vanishing (static) octupole mo- 
ments (see, Tables I and II). Nevertheless, after both projection onto it = +1 and dynamical 
(Q20; <33o)-fluctuations are considered at the 2D-GCM level, the octupole deformation effects 



24 



predicted for 148 > 150 Sm and 150 > 152 Gd are reduced to more than half of their mean field values. 
At variance with the HFB results, the nuclei 146 ' 152 > 154 Sm and 148 > 154 . 156 Gd exhibit dynamical 
ground state octupole moments Q$ ~ 0.40-0.50 b 3 / 2 . We conclude that, regardless of the 
particular version of the Gogny-EDF employed, our 2D-GCM calculations suggest a dynam- 
ical shape/phase transition from weakly ( 146 Sm and 148 Gd) to well quadrupole deformed 
( 154 Sm and 156 Gd) ground states as well as a transition to an octupole vibrational regime in 
the considered Sm and Gd nuclei. 

For the lowest-lying negative parity states, the dynamical octupole and quadrupole 
Q20 moments, computed with the corresponding 2D-GCM states l^Z^ 1 } or 1^=3 1 ), are 
also listed in Tables III and IV. It should be noted that the largest values of the octupole 
deformations Q$ and always correspond to the N=88 isotones 150 Sm and 152 Gd. 

The values of the ground state dipole moments are less predictable than the averages 
of the quadrupole and octupole moments discussed previously as the behavior of D for the 
HFB states strongly depends upon the orbitals occupied and those change rapidly with 
deformation. The comparison of the dipole moments with available experimental data [l| 
is presented in panels a) and b) of Fig. [7J In particular, the comparison between the HFB 
results (see, Tables I and II), Dq + ^ and experimental values clearly reveal the limitations of 
the HFB approximation to predict dipole moments in this region of the nuclear chart. 

Another physical observable is the energy splitting between the lowest lying 7r = +1 
and Ti = — 1 states. The results for 146-154 Sm and 148_156 Gd are compared in Fig. [8] with 
available experimental + — 1~ and + — 3~ energy splittings [74]. As already mentioned, in 
the present study we are not able to account for genuine 1~ and/or 3 _ states that require, for 



example, the use of cranking HFB states 3l|, |4l| . With this in mind and regardless of the 



Gogny-EDF employed, a reasonable agreement between the theoretical and experimental 
energy splittings is observed. The remaining discrepancies imply that correlations other 
than (axial) quadrupole-octupole fluctuations could also be required. In particular, the 
time-odd components of the Gogny-EDF, incorporated throughout cranking calculations, 
should be further investigated within the present 2D-GCM framework. Let us mention 



that the results are compatible with the ones obtained in Ref. 



33] using a one-dimensional 



collective Hamiltonian whose parameters are derived from octupole constrained calculations. 



This is also the case with the systematic calculations of Ref. 
octupole degree of freedom as generating coordinate. 



60| using a GCM with the 
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In panels a) and b) of Fig. [9j the B(E1, 1 — > + ) reduced transition probabilities of Eq. 
IT]) are compared with experimental data [1]. It is very satisfying to observe how, without 



resorting to any effective charges, the predicted B(E1, 1~ — > + ) values in Sm nuclei follow 
the experimental isotopic trend with a slight improvement in the case of the Gogny-DIM 
EDF. In panels c) and d) of the same figure, we compare the B(E3,3~ — > + ) transition 



rates of Eq. (JT7J] with available data [73j]. The predicted B(E3, 3~ — > + ) values reproduce 
quite well the experimental ones in the case of 152,154 Sm and 154,156 Gd. On the other hand, 
from the comparison between ours and the B(E1,1~~ — > + ) and B(E3,3~ — > + ) rates 



obtained in Refs. |33j and |60j, we can conclude that they are, to a large extent, not very 
sensitive to quadrupole fluctuations. 

In panels a) and b) of Fig. (TUl the correlation energies defined as the difference between 
the reference HFB ground state energy and the 2D-GCM one 

^corr - ^HFB ~ Ar=l l 2U J 

are ploted. The parity restoration correlation energies E^orr of Eq. ([8]) are also included 
for comparison. The predicted isotopic trends and quantitative values of E^®~ GCM are quite 
similar for both Gogny-DlS and Gogny-DIM EDFs. The correlation energies E%®~ 
exhibit a relatively weak dependence with neutron number with values oscillating between 
1.74 and 2.09 MeV for Sm and between 1.83 and 2.17 MeV for Gd nuclei. The smooth 
variation of the correlation energy is, however, of the same order of magnitude as the rms 
for the binding energy in modern nuclear mass tables [4^ and therefore the dynamical 
octupole correlation energies should be considered in improved versions of them. 

A rough estimate of the contribution of the (Q2Cb Q30) fluctuations to the correlation 
energies can be obtained by subtracting to the total correlation energy the parity projected 
one. Those contributions range in between 0.94 and 1.41 MeV for Sm isotopes and in 
between 0.94 and 1.56 MeV for Gd isotopes. The oscillations are slightly larger than for the 
total correlation energy. 

In order to determine the contributions of each degree of freedom in the results obtained 
we have also performed one dimensional GCM calculations along each of the degrees of 
freedom. First, the octupole moment has been used as generating coordinate. For each 
octupole moment considered the quadrupole moment corresponds to the minimum energy. 
The octupole moments of the generating wave functions are taken in the range -7 b 3 / 2 
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< Qso < 7 b 3 / 2 and with a mesh size SQ 30 = 0.25 b 3 / 2 . The 1D-GCM ansatz is 

\Kid-qs) = J dQnftiD-wiQnmQx)) (21) 

given in terms of the HFB states |$(Q3o))- Note that no quadrupole constraint is imposed in 
these calculations. From the 1D-Q3 ground state energies we can computed the ID octupole 
correlation energy 

^corr — J- 1 HFB ^cr=l,lD~Q3 \ ZZ ) 

This quantity is displayed in panels a) to d) of Fig. [11] for the considered Sm and Gd nuclei. 
It has to be mentioned that this type of calculations have been carried out for all possible 



even-even nuclei with several parametrizations of the Gogny force in [60]. 

In a second step, GCM calculations with the quadrupole degree of freedom (Q 30 = 0, 
i.e., reflection symmetry is preserved) as generating coordinate have been performed. The 
Q20 values used are in the interval -30b < Q20 < 30b with 6Q 2 o= 0.6 b. The GCM wave 
functions 

\**,1D-Q2) = J dQ 20 fa,lD-Q2(Q20)mQ20)) (23) 

are defined in terms of the states \$(Q2o))- The corresponding correlation energy 

E lorr Q2 — E HFB ~ E(T=1,1D-Q2 (24) 

is displayed in panels a) to d) of Fig. [TU 

In panels a) to d) of Fig. HH we compare the sum E^ Q2+1D - Q3 = E^- Q3 + E\®~ Q2 
with the correlation energies E 2 ® r ~ GCM of Eq. (|20|) . For the particular set of Sm and 
Gd nuclei considered in the present study and regardless of the Gogny-EDF employed, the 
correlation energies provided by the full 2D-GCM are very well reproduced by the sum of the 
ones obtained in the framework of the 1D-GCM approximations (j2Ti) and ( |23l) . Obviously, 
this is far from being a general statement and further explorations in other regions of the 
nuclear chart, specially those showing shape coexistence already at Q30 = are required. 
Nevertheless, the kind of decoupling observed in our results may be potentially relevant to 
incorporate correlation energies stemming from parity restoration and octupole fluctuations 
in large scale calculations of nuclear masses based on the Gogny-EDF (see, for exam ple, 



53l |) as well as in future fitting protocols beyond the most recent DIM parametrization 



43] 



of the Gogny EDF. 
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V. CONCLUSIONS 



Calculations have been carried out using the GCM method and with the multipole mo- 
ments Q20 and Q 30 as generating coordinates for several Sm and Gd isotopes and with dif- 
ferent parametrizations of the Gogny force. The results from different parametrizations are 
very close to each other indicating again that the DIM parametrization of the Gogny force 
performs as well as D1S in spectroscopic calculation. The comparison with experimental 
data is fairly good both for excitation energies and electromagnetic transition probabilities 
reassuring the predictive power of the Gogny class of EDFs. Comparison of the 2D GCM 
results with the outcome of previous ID Collective Schrodinger equation calculations in the 
same region points out to a decoupling of the dynamics of the quadrupole and octupole 
degrees of freedom. This conclusion is reinforced by the comparison of the 2D correlation 
energies with the sum of correlation energies along each of the degrees of freedom. Cor- 
relation energies show a smooth behavior with neutron number with differences between 
different isotopes as large as 200 keV. Although these differences are small, they can be rele- 
vant for theories aiming at providing accurate mass tables for applications requiring accurate 
reaction rates that depend on their energetic balance. 
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